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ABSTRACT 

The Rossiter-McLaughlin (RM) effect is the distortion of stellar spectral lines that occurs during 
eclipses or transits, due to stellar rotation. We assess the future prospects for using the RM effect 
to measure the alignment of planetary orbits with the spin axes of their parent stars, and to confirm 
exoplanetary transits. We compute the achievable accuracy for the parameters of interest, in general 
and for the 5 known cases of transiting exoplanets with bright host stars. We determine the require- 
ments for detecting the effects of differential rotation. For transiting planets with small masses or 
long periods (as will be detected by forthcoming satellite missions), the velocity anomaly produced 
by the RM effect can be much larger than the orbital velocity of the star. For a terrestrial planet 
in the habitable zone of a Sun-like star found by the Kepler mission, it will be difficult to use the 
RM effect to confirm transits with current instruments, but it still may be easier than measuring the 
spectroscopic orbit. 

Subject headings: planetary systems — stars: rotation 



1. INTRODUCTION 

When an exoplanet transits the disk of its parent star, there is both a photometric signal and a spectroscopic signal. 
The photometric signal is a small reduction in the received flux due to the partial obscuration of the stellar disk, 
as first detected by Charbonneau et al. (2000) and Henry et al. (2000). The transit light curve depends chiefly on 
the radius of the planet, the radius of the star, the orbital inclination, and the stellar limb darkening function. The 
spectroscopic signal is subtler and less familiar. Changes in the stellar absorption lines can be produced by absorption 
features in the planetary atmosphere (see, e.g., Charbonneau et al. 2002, Vidal-Madjar et al. 2003). Or, if the emergent 
stellar spectrum varies with position across the stellar disk, then the planetary obstruction will cause changes in the 
disk-integrated spectrum. 

In particular, the Rossiter-McLaughlin (RM) effect refers to the spectral distortion caused by the spatial variation in 
the emergent spectrum due to stellar rotation. The exposed portion of the photosphere has a net rotational Doppler 
shift. When the planet covers part of the blueshifted half of the stellar disk, the integrated starlight appears slightly 
redshifted, and vice versa. Thus, the spectral distortion of the RM effect is often manifested as an "anomalous" radial 
velocity, i.e., a Doppler shift that is greater or smaller than the shift expected from only the star's orbital motion. 
Although the RM effect has a long history in the context of eclipsing binary stars (Forbes 1911, Schlesinger 1911, 
Rossiter 1924, McLaughlin 1924), its importance in the context of exoplanets is only starting to be appreciated. 

Suppose that the photometric signal has already been observed, and that the parameters governing the photometric 
signal have been accurately determined. The two most important additional parameters that govern the RM effect are 
Vs sin/g, the projected rotation speed of the stellar surface, and A, the angle between the sky projections of the stellar 
spin axis and the orbit normal (i.e., the angle between the transit chord and lines of stellar latitude, for Is = 90°). 
Observations of the RM effect allow these two parameters to be measured (see, e.g., Queloz et al. 2000, Ohta et al. 
2005, Winn et al. 2005, Wolf et al. 2006). The latter parameter, A, is especially interesting, as it provides information 
about exoplanetary spin-orbit alignment. Solar system planetary orbits are generally within ~5° of the solar equatorial 
plane (Beck & Giles 2005), but this may or may not be true for the full range of exoplanetary systems. 

Of particular interest in this regard are hot Jupiters — planets with masses of ~A/j up and orbital periods of ^3 days — 
because there are some theoretical reasons why one might expect large misalignments for such planets. Hot Jupiters are 
thought to have formed at large orbital distances, and then migrated inward to their current positions. The migration 
mechanism is still unknown, and some of the proposed mechanisms differ in the degree to which they would affect 
spin-orbit alignment. Thus, measurement of spin-orbit alignment offers a possible means of discriminating among 
migration theories. The most widely discussed category of migration theories involves disk-planet interactions (Lin 
et al. 1996). There is a large literature on this subject (as reviewed by Papaloizou & Terquem 2006), but broadly 
speaking, there is no apparent reason why these mechanisms would perturb spin-orbit alignment; in fact they may 
even drive the system toward closer alignment. In contrast, other migration theories involve disruptive events such 
as planet-planet interactions (Rasio & Ford 1996, Weidenschilling & Marzari 1996) or planetesimal collisions (Murray 
et al. 1998), which could act to randomize spin-orbit alignment. Another proposal involves the Kozai mechanism, in 
which a companion star causes oscillations in the planetary orbit's eccentricity and inclination (Innanen et al. 1997, 
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Holman et al. 1997). By the time tides circularize the orbit and halt "Kozai migration," the orbital inclination can 
change substantially (Wu & Murray 2003; Eggenberger, Udry, & Mayor 2004; D. Fabrycky & S. Tremaine, priv. 
comm.). 

For planets at larger orbital distances, there is no particular reason to expect large misalignments, but the field of 
cxoplanets has rewarded observers with surprises in the past. There are also a few empirical hints of multiple orbital 
planes in some systems, such as the double debris disk recently reported around (3 Pic (Golimowski et al. 2006), and 
the apparently counter-rotating disks around the protostar IRAS 16293-2422 (Remijan & Hollis 2006). 

The exoplanetary RM effect was first observed by Queloz et al. (2000) and Bundy & Marcy (2000) during transits 
of HD 209458b. The former authors were able to place an upper bound on A of about 20°. Snellen (2004) described 
the interesting idea of using the wavelength-dependence of the RM effect to search for planetary absorption features, 
but with only null results for HD 209458b. Winn et al. (2005) used improved photometric and spectro scopic data 
to sho w that A = — 4?4 ± 1?4, a small but significant misalignment reminiscent of Solar system planets. iWolf et all 
(2006) performed a similar study of HD 149026b, finding A = 11 ± 14°. Apparently, in these two cases, the migration 
mechanism was fairly quiescent. We can soon expect similar stud ies of the other tran siting pla nets with host star s 
bright enough for high- precision spec troscopy, namely, HD 189733b ijBouchv et all2005cj) . TrES-1 l|Alonso et all2004[) . 
and XO-lb (McCullou gh et aLl l2006h The first goal of this paper is to provide useful guidance for future observations 
of these and other systems. 

The second goal of this paper is to investigate another important and timely application of the RM effect: transit 
detection and confirmation. T he spectroscopic dete ction of transits offers certain advantages, in some cases, over 
photometric detection. Indeed, Bou chv et alJ ( 2005^ discovered the transits of HD 189733b through observations of 
the RM effect, and followed up this discovery with photometric observations. The reason for this ordering may simply 
have been that the observers had more convenient access to spectroscopic resources than to photometric resources. 
However, as pointed out by Welsh et al. (2004) and Ohta et al. ( 2005), the im portance of the RM effect in tr ansit 
confirmation is likely to increase soon. The satellite missions Corot (Baglin 2003) and Kepler ijBorucki et alJl2003|) aim 
to find smaller and longer-period transiting planets than those currently known, for which it will be very difficult to 
confirm the existence of transits photometrically from the ground. It has been envisioned that the confirmation process 
will require measuring the orbital velocity of the parent star, and there is ho pe that this can be achieved, given th e 
recent excellent progress in improving the accuracy of Doppler measurements l)Mavor et al.ll2003 iMarcv et alj [2005) . 
As we will show in this paper, the RM effect offers an alternative path to confirmation that is suitable for at least a 
subset of stars. 

Previous theoretical work on this topic has concentrated on analytic descriptions of the spectral distortion (Ohta 
et al. 2005, Gimenez 2006), or on numerical simulations of the spectral distortion involving a discretized model of 
the stellar surface (Queloz et al. 2000, Welsh et al. 2004). In this work, we are not concerned with high accuracy 
in describing the spectroscopic distortion; instead we are concerned with the measurement and estimation problem. 
In § |5J we explain our notation, review some of the previously derived results and provide some useful approximate 
scaling relations. In §0 we estimate the achievable accuracy in measuring the key RM parameters, Vs sin Ig and A, as 
a function of the orbital geometry of the system and of the characteristics of the data. We do this both for the general 
case, and for the specific cases of the 5 transiting exoplanets with bright host stars. A secondary parameter that 
affects the RM signal is the degree of differential rotation across the stellar disk, and in §0]we investigate whether or 
not this effect is important for near-term observations. In § we derive an analytic expression for the signal-to-noise 
ratio (S/N) in the detection of the RM effect, and we apply this formula to assess the prospects for state-of-the-art 
Doppler measurements to confirm the existence of transits of planets detected by Kepler and Corot. The final section 
summarizes all of these results and suggests some avenues for future work. 

2. THE ROSSITER-MCLAUGHLIN EFFECT 

A pedagogic illustration of the physics of the Rossiter-McLaughlin effect is presented in Figure ^ The top row of 
panels illustrates the progress of a transiting planet across the limb-darkened disk of its parent star. This is the origin 
of the photometric signal. In the second row of panels, the stellar disk has been color-coded according to the local 
line-of-sight velocity due to rotation. One half is blueshifted, and the other half is redshifted. This is the origin of the 
RM signal. The third row of panels is a schematic view of a stellar absorption line that would be recorded in unresolved 
(i.e., disk-integrated) observations, for the case in which rotation is the dominant source of line broadening. The planet 
hides a small fraction of a small range of velocity components. The result is a "bump" that moves through spectral 
line as the planet moves across the star. The idea is essentially the same one that is employed in Doppler imaging 
of rapidly rotating stars (see, e.g., Rice 2002), except in this case the contrast is provided by an opaque foreground 
object rather than star spots. The final row of panels depicts the case in which rotation is not the dominant source of 
broadening. The distortion in this case is most easily detected as a net Doppler shift. 

Most of the known exoplanets have been discovered by searching for periodic Doppler shifts in the spectrum of a 
star that result from the star's reflex velocity (i.e., its orbital motion around the star-planet center of mass). The 
target stars in these searches are deliberately chosen to be inactive main-sequence FGK stars, which by nature are 
slow rotators ( ,$ 5 km s"" 1 ). The reason is that those stars offer the smallest intrinsic velocity noise. This means that 
the RM effect is appropriately described as an anomalous radial velocity. When a spectrum taken during transit is 
compared with a standard stellar template in order to measure the radial velocity shift, the distortion of the spectral 
lines will induce a signal which appears as an anomalous radial velocity. 

How can the size of the anomalous velocity be predicted, in terms of the planetary and stellar properties? If, for 
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Fig. 1. — The physics of the RM effect. Top row. Three successive phases of an exoplanetary transit. Second row. Same, but 
the projected stellar rotation speed at each point has been color coded, neglecting differential rotation. At each phase, the planet hides a 
different velocity component. Third row. Illustration of an observed stellar absorption line, for the case of purely rotational broadening, 
i.e., the net broadening due to all other mechanisms is much less than the rotational broadening (Wp -C V^sin/g). The missing velocity 
component is manifested as a time- variable bump in the line profile. Fourth row. Same, but for the case W p ~ Vgsin/g, in which other 
line-broadening mechanisms besides rotation are important. 



example, the observed spectrum is simply cross-correlated with a standard stellar template, then there would be a 
wavelength shift in the peak of the cross-correlation function relative to the unocculted case. In this case, an excellent 
approximation for the velocity anomaly can be obtained by computing the first moment of the distorted line profile 
(i.e., the flux- weighted mean wavelength), and comparing it to the first moment of the line profile in the out-of-transit 
spectrum, as was done by Ohta et al. ( 2 005) and Gimenez (2006). 

In fact, as discussed bv I Winn et al.l {2005), the accuracy of this approximation depen ds on the exact pro cedure 
by which radial velocities are extracted from a series of observed spectra. For example, iButler et alJ lfl99 61 use a 
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procedure that is optimized for high precision velocity measurements using an iodine reference cell. It is considerably 
more complicated than a simple cross-correlation, and has as its working assumption that all spectral cha nges are due 
only t o an overall Doppler shift (which is to be measured) and variations in the focus or the instrument. I Winn et alJ 
(2005) gave evidence that in this case the first-moment approximation for the RM effect is only accurate to within 
~10% for HD 209458. This level of accuracy is sufficient for our main purpose of providing scaling relations and 
estimated measurement accuracies within a factor of 2. For this reason, we will employ the first-moment approximation 
throughout this paper. However, it is interesting to note that the reported size of the anomaly will depend on the 
precise method by which radial velocities are extracted, which may make it difficult to compare the results of different 
investigators even when they observe the same system. 

Consider a star with mass M and radius R that has a transiting planet with mass m and radius r. The orbit has 
a period P, eccentricity e, and argument of pericenter u>. We can write the net radial velocity variation of the star 
{AV) as the sum of the radial velocity of the star due to its orbital motion (AVb) and the anomalous radial velocity 
due to the RM effect (AV R ): 

AV(t) = AV (t) + AV R (t). (1) 

The "O" stands for orbit, and the "R" refers to both Rossiter-McLaughlin and rotation. 
The line-of-sight component of the orbital velocity is 

AV (t) = K {cos[f(t) +uj] +ecostj}, (2) 

where / is the true anomaly and Ko is the orbital velocity semiamplitude, 

Ko= {—) (M + w)2 /3 ( X - e ) ( 3 ) 
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where / is the orbital inclination with respect to the sky plane. In the latter equality we have assumed m -C M and 
e = 0. 

Assuming that the width of the absorption line is dominated by rotational broadening, and further assuming that 
the stellar Doppler shift is small, the first-moment approximation mentioned previously gives (|Ohta et alJ l2005): 

AV R (t) - - V s sml s — . (4) 

J JZ{x, y)dxdy 

Here, Vs is the equatorial rotation speed of the stellar photosphere, Is is the inclination of the stellar spin axis relative 
to the sky plane, and T(x, y) is the surface brightness of the observed stellar disk (including the dark spot due to 
the planet). The sky-plane coordinates x and y are measured in units of the stellar radius, have their origin at the 
projected center of the star, and are perpendicular and parallel to the projected stellar rotation axis, respectively. 
In fact, equation also holds for lines that have additional broadening mechanisms, such as thermal broadening, 
provided that the additional broadening mechanisms produce no net Doppler shift (i.e., the broadening kernel is 
symmetric about its centroid). 
For convenience, we write the RM effect as 

AV R (t) = K R g(t;x P) y p) j,e,...), (5) 

separating the overall amplitude K R of the RM effect from the dimensionless function g(t) < 1. The amplitude is 
given by 

7 2 
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K R = V s sin I s J- , (6) 
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where 7 = r/R. In the latter equality, we have assumed 7 <C 1. For convenience, we will define V = Vs sin Is- The 
dimensionless function g depends primarily on the projected position of the planet (x p ,y p ), but also on 7 and the 
limb-darkening function. For simplicity, we use a single-parameter "linear" description of the limb darkening law, such 
that the (unocculted) surface brightness of the star is 



1{x,y) _ 
2o 



l-(l-.* 2 -y 2 ) 1/2 



(7) 



with e the linear limb darkening parameter. Note that in some circumstances — for example, the case of differential 
rotation as discussed in ^ — the function g will depend on additional parameters. 

Figure |21 shows three different trajectories of a transiting planet across the stellar disk. These trajectories all have 
the same impact parameter b, and consequently they all produce exactly the same photometric signal. 3 However, the 

3 The impact parameter is given by b = a cos I /R, where a is the orbital semimajor axis. 
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Fig. 2. — The dependence of the RM waveform on A. Three different possible trajectories of a transiting planet are shown, along with 
the corresponding RM waveform (as computed with the formulae of Ohta et al. 2005). The trajectories all have the same impact parameter 
and produce the same light curve, but they differ in A and produce different RM curves. The dotted lines are for the case of no limb 
darkening (e = 0), and the solid lines are for e = 0.6. 



trajectories differ in the value of A, and consequently produce different RM waveforms, as plotted in the lower row 
of panels. The sensitivity of the RM waveform to A is what enables the observer to assess spin-orbit alignment. The 
question of the achievable accuracy in A will be taken up in the next section. 

An especially simple case is when the planetary disk is fully contained within the stellar disk, and limb darkening 
is negligible (e = 0). In that case, g is the perpendicular distance from the projected stellar spin axis, g(t) = x p (t). If 
we consider a rectilinear trajectory across the face of the star with impact parameter 6, we can write the position of 
the center of the planet as a function of time as, 

x p (t) — t cos A — b sin A, 

y p [t) =rsin A + 6 cos A, (8) 

where t = (t—tt Ia )/Tt Ia , t tra is the time of the transit midpoint, T tla — i?/u or b is the radius crossing time corresponding 
to the planet's orbital velocity at the time of transit (so that the transit duration is approximately 2Tt ra yl — W) and 
A is the angle of the trajectory with respect to the the apparent stellar equator. We define A to be between —180° and 
+ 180°, such that for A > 0, the planet moves towards the stellar north pole as it proceeds across the stellar disk. We 
note that due to the rotation of the star, the familiar symmetry b «-> — b of the photometric signal is broken, and thus 
it is important to specify the sign of b when considering the RM effect. 

From equation ©, it is clear that in the case of uniform rotation, no limb-darkening, and when the planet is fully 
contained within the stellar disk, the RM curve is simply a linear function of time. During ingress and egress, and 
when limb darkening is taken into accou nt, t he expres s ion fo r g is more complicated and may need to be evaluated 
numerically. However. lOhta et alJ (|2Q05) and lCimenea ( 2006) provide useful approximate analytic expressions. 

For our order-of-magnitude calculations, we will assume that the planet is small (7 <C 1), and we will not be 
concerned with the ingress and egress phases because they constitute only a small fraction (~27) of the entire duration 
of the transit. For now, we will also neglect limb darkening. We therefore adopt the simple approximation 

a(t) = { X p W if Vt 2 + b 2 < 1, fQ) 
yw 1 otherwise. v ' 

It is interesting to compare the amplitude of the stellar orbital velocity Ko with the amplitude of the anomalous 
velocity Kr. Assuming circular orbits (e = 0), 7 <C 1, and / = 90°, we find 

K \2TTGm) \pj ' {W) 

where p p and p* are the average densities of the planet and star, respectively. Since we expect p*/p p to vary by no 
more than a factor of a few from system to system, the order of magnitude of Kr/Ko depends mainly on the orbital 
period, the mass of the planet, and the projected rotation speed of the star. We find 

^-03 (^-] V3 ( V ) (11) 

Ko ' \M Jup J Vsdays/ \5 kms' 1 J V 1 
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All of the currently-known transiting exoplanets have masses that are comparable to Jupiter's mass, and orbital 
periods of 1-4 days. For these systems, the anomalous velocity is smaller than the orbital velocity by a factor of a 
few. However, the properties of the known systems have been subject to very strong selection effects: the transits of 
large, short-period planets are much easier to detect than those of small, long-period planets (Gaudi 2005, Gaudi et 
al. 2005). An interesting implication of equation (jl 1|) is that for the most challenging systems (small planets, long 
periods), the amplitude of the RM effect will exceed the stellar orbital velocity. In particular, for an Earth- mass planet 
with a period of one year, Kr/Ko ~ 3 for V = 5 kms -1 . This explains the appealing possibility of using the RM 
effect to detect or confirm transits, which will be discussed further in § 

3. PROSPECTS FOR MEASURING SPIN-ORBIT ALIGNMENT 

In this section, we consider what can be learned from observations of the RM effect in the regime of large S/N. We 
have in mind high-cadence, high-precision observations of both the photometric and spectroscopic transit by a hot 
Jupiter, for which the spectroscopic orbit is already well established. 

As mentioned in § 1, a principal goal of such studies is the determination of A, the angle on the sky between the 
stellar spin axis and the planetary orbit normal, because this angle gives a lower bound on any misalignment between 
the angular momenta of the star and the orbit. The value of A is determined or bounded by fitting a parameterized 
model to the photometric and spectroscopic measurements (see, e.g., Winn et al. 2005 and Wolf et al. 2006). However, 
for simple calculations, and for planning purposes, it is useful to have some heuristics and order-of-magnitude estimates 
for the effect of A on the RM waveform. 

In particular, it is useful to consider the relative timing of three observable events: the moment of greatest transit 
depth (ttra); the moment when the orbital radial velocity variation is zero (t rb)j and the moment when the anomalous 
RM velocity variation is zero (t ro t)- At t = tt ra , the projected planet-star distance is smallest. At t = i rbj the star is 
moving in the plane of the sky. At t = t rot , the planet lies directly in front of the stellar rotation axis. For a circular 
orbit with A = 0, these three events are simultaneous. If the orbit is circular but A ^ 0, then t tra = Wb but t Iot will 
occur either earlier or later: 

Ttra 

This is easily derived from equation JSJ). For a noncircular orbit, ttra and t OT b are no longer simultaneous. To first 
order in e, 

P 

torb ~ ttra = 7^ ( e C0S W ) I ( 13 ) 
Z7T 

where uj is the argument of pericenter. To the same order, both the time difference i ro t — itra and the transit duration 
Tt ra are multiplied by the same factor 1 — esinw, and hence Ea. 1121 remains valid. These expressions give some sense 
of the timing accuracy that is needed for a desired accuracy in A. For example, for a mid-latitude transit at b = 0.5 
lasting 2.5 hr, a misalignment of A ~ 1° corresponds to a timing offset of 45 seconds between the transit midpoint and 
the null in the RM waveform. 

Next, we derive an expression for the expected uncertainty in A based on a series of spectroscopic measurements. 
Consider a series of N radial velocity measurements taken at times tk during the planetary transit (between first and 
fourth contact), each of which has an uncertainty a. We consider the case in which the RM measurements are the 
limiting source of error; we assume that both the photometric transit and the spectroscopic orbit have already been 
measured accurately. Thus the times of contact, t t ra, and the parameters 7 and b have negligible uncertainties, and 
the orbital velocity AVo(t) can be accurately subtracted from the total velocity variation observed during transits to 
isolate the RM waveform. The only parameters to be determined by fitting the transit data are A and V, which we 
combine into a two-dimensional parameter vector a. The expected uncertainties in A and V are the square roots of 
the diagonal elements of a matrix C that is given by 



C = B- 



where B is related to the Fisher information matrix fsee IGouIdll2003|) and is calculated as 



N 1 



fc=i 



dai 



^-AV R (t k ] 



(14) 



(15) 



We will adopt the approximate analytic form for AV R (t) that was given in equation @. Assuming evenly-spaced 
observations and large N, we can convert the sum in equation (|15fl into an integral. This then yields the following 
expressions for the expected uncertainties in A and V: 

"(l-& 2 )sin 2 A + 3& 2 cos 2 A n 1/2 
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where wc have defined 
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(16) 
(17) 

(18) 
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As discussed further in SjSJ the factor Qr is proportional to the total signal-to-noise ratio of the measured RM 
waveform. Note that <jy(X)/V = <j\(tt — A). Similarly, we can derive the covariance between A and V, which is given 
by cov(X,V) = Ci 2 /(CiiC 2 2) 1/2 . We find 

46 2 - 1 

cov(A, V) = m . (19) 

[(1 - 6 2 ) 2 + % 4 + 36 2 (1 - & 2 )(tan 2 A + cot 2 A)] 1 

The uncertainty in each parameter is the product of Q^ 1 and a factor that depends on the orbital geometry. The 
geometrical factor is illustrated in Figure |21 for 7 = 0.1 and three values of the impact parameter b. The solid lines 
show the expected uncertainties in A and V given by equations l|16[) and l|17|l. after dividing by QZ . The three cases 
are a near-central transit (b = 0.1), a mid-level transit (b — 0.5), and a grazing transit (b = 0.9). For b = 0.5, we 
find that o\ — ay /V = 2Q~j, , and that the uncertainties in V and A are uncorrelated [cov(X,V) = 0]. In this sense, 
mid-level transits are ideal for cleanly separating the effects of Vs sin/5 and A on the RM waveform. 

In contrast, for nearly central transits there is a strong degeneracy between A and V and a strong dependence of a\ 
on A. For b = exactly, cov(A, V) = — 1. The observed signal depends only on the parameter combination V cos A, and 
the timing offset given in Eq. (|13fl is zero regardless of A. One can measure only the amplitude of the RM waveform, 
and cannot tell whether a small amplitude (say) is the result of an equatorial transit across a slowly rotating star, or 
a misaligned transit across a more rapid rotator. An accurate measurement and interpretation of the line broadening 
in the out-of-transit spectra is essential here, by providing an independent estimate of V that can be used as an a 
priori constraint. For grazing transits, there is a more modest dependence of the uncertainties on A (although our 
approximations are least accurate for grazing transits, as described below). 

In order to verify our analytic expressions, and evaluate the importance of some of the effects that these expressions 
neglect (namely, the finite ingress/egress durations and limb darkening), we numerically eval uated the elements of the 
Fisher matrix B using the more accurate but more complex expressions for AVr(£) given bv lOhta et alJ ((2005). The 
results are also plotted in Figure |3 The dotted lines show the results for e = (no limb darkening), and the dashed 
lines show the results for e = 0.6. Except for non-grazing transits, the differences between the uncertainties predicted 
by our analytic expressions and the numerically-calculated uncertainties are small (< 10%). Not surprisingly, for 
grazing transits the differences are substantially larger and can be a factor of 2. 

Next, we verified our results using Monte Carlo simulations. For a given choice of V and A, we created 5000 simulated 
data sets of radial velocity measurements during a transit. Each data point was the value calculated according to the 
expressions of Ohta et al. (2005), plus Gaussian noise with a standard deviation of a. We then used a downhill-simplex 
algorithm to optimize the parameter values V and A for each data set. We calculated the dispersions in the resulting 
distribution of 5000 fitted values, and took the dispersion to be the uncertainty in that parameter. We repeated this 
procedure for a range of V and A. The resulting uncertainties agree very well with the uncertainties computed by 
numerically evaluating the elements of B as described above. 

We now consider the 5 particular cases of known transiting exoplanets with bright ( V < 12) parent stars : HD 209458 
(Henry et al. 2000, Charbonneau et al. 2000), TrES-1 (Alonso et al. 2004), HD 189733 l)Bouchv et all2005cft . HD 149026 
(Sato et al. 2005), and XO-1 (McCullough et al. 2006). For the systems for which the RM effect has already been 
measured, this exercise provides an empirical check on our calculations. For the others, it provides a guide for the 
achievable accuracy of future observations. Below we describe our assumptions and the results for each planet in 
detail; this information is summarized in Tabled 

For HP 209458, an analysis of the best available photometric and spectroscopic data was performed bv lWinn et alJ 
(2005). The data set included N = 14 radial velocity measurements during transits, with a typical uncertainty of 
a = 4.1 ms . By fitting a parameterized model to all of the data, and evaluating the parameter uncertainties with a 
Monte Carlo bootstrap technique, they found A = — 4?4± 1?4 and V = 4.70±0.16 kins" 1 . Using these values, as well 
as the parameters 7 = 0.121 and b = 0.52 determined from the photometry, our analytic estimates predict that the 
uncertainties should be o\ = 177 and ay = 0.15 kms -1 . These estimates are in excellent agreement with the actual 
uncertainties derived from detailed model-fitting. In addition, Eq. Ijl9|) predicts that the uncertainties in V and A 
should be uncorrel ated, as was i ndeed found to the be the case by Queloz et al. (2000) and Winn et al. (2005). 

For HD 149026, iWolf et all l)2006[) analyzed all of the available photometric and spectroscopic data, including 
N = 15 radial velocity measurements taken during transit. The typical velocity uncertainty was a — 4.0 ms . 
Through parameterized model-fitting and bootstrap resampling, these authors found A = 11 it 14° V = 6.4^ 'y kms -1 , 
7 = 0.0 51, and b = 0.39. Our analytic formulae predict a\ = 9.0°and ay = 0.76 kms -1 . Our estimates are smaller 
than the Wolf et alJ l)2006|) uncertainties. We believe that the reason is that the photometric signal is not nearly as well 
determined for HD 149026 as it is for HD 209458, not only because of the smaller size of the planet, but also because 
the existing observations are from ground-based telescopes as opposed to the higher-precision Hubble Space Telescope 
light curve obtained by Brown et al. (2001) for HD 209458. Consequently, the uncertainties in the parameters ftra,7, 
and b cannot be neglected. The covariance with these photometric parameters contribute to the uncertainty in A and 
V, violating one of the conditions for the accuracy of our analytic estimates. We conclude that there is scope for 
improvement in the determination of A for this system through improved photometry. 

As mentioned in the introduction, the transits of HD 189733b were originally discovered via the RM effect 
(jBouchv et alJl2005cj) . but those authors did not attempt to measure A or V by fitting the RM waveform. However, 
using their estimate of V = 3.5 ± 1.0 kms -1 derived from the out-of-transit spectral line profile, we can anticipate the 
expected uncertainties of such attempts, through an application of our expressions. The photon noise in the radial 
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TABLE 1 

Predictions for. Rossiter-McLaughlin Effect Parameter Uncertainties 
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Note. — For HD 189733b, a fit to the Rossiter-McLaughlin effect data has not been reported, and so no constraint on A is available. 
We therefore show the expected parameter uncertainties for the two values A = 0, 90° , which will bracket the range of uncertainties. For 
TrES-lb and XO-lb, no Rossiter-McLaughlin data has been reported, and so we assume N = 15, a = 15 ms -1 , and show the expected 
parameter uncertainties for the two values A = 0, 90° . 

a V = V s sin/ S 

b For TrES-lb, we assume a projected stellar velocity that is equal to the spectroscopically determined upper limit of V = 5 kms -1 . 



velocity measurements reported bv IBouchv et al.l l|2005d) amounts to only ~5-7 ms -1 , but the actual uncertainties 
are larger because the star is fairly active. As an estimate of this "stellar jitter," we will adopt the value a = 15 m s" 1 
based on the scatter around the Keplerian orbital solution. Taking b = 0.66 and 7 = 0.157 (B akos et al.ll2006l) . we 
predict a\ — 5?2 and ay = 0.48 kms -1 for N = 8 (again, assuming that the measurement of the RM waveform is the 
limiting uncertainty). These values are based on the assumption A ~ 0, although the expected uncertainties are fairly 
insensitive to A because the transit is neither central nor grazing (see the previous section). If the velocity jitter turns 
out to be smaller, or if a larger number of transit velocities are measured, then the uncertainties in the RM parameters 
will scale as a j\J~N . 

No measurements of the RM effect have yet been reported for TrES-1 or XO-1. We estimate the accuracy with which 
A can be constrained in these systems with future observations, using the parameters of the systems that have been 
derived from the existing photometric and spectroscopic data. For specificity, we assume N = 15 and a = 15 ms . 
The latter estimate for a seems reasonable, given that the apparent magnitudes of the host stars (my = 11 — 12) are 
fainter than the three stars considered above (my = 8). 

For TrES-1, there exists only an upper limit on the projected stellar rotation velocity of V < 5 kms -1 . We will 
optimistically assume that the true value is 5 kms -1 , right at the upper limit. Then, assuming b = 0.18 and 7 = 0.137 
ijAlonso et al.ll2004t iSozzetti et aTll2004j) . we find a x = 13° and cry = 0.36 kms -1 for A = 0, and cr A = 4?1 and 
ay = 1.1 kms -1 for A = 90°. The strong sensitivity to A arises because the transit is nearl y central. 

For XO-1, we adopt V = 1.11 kms -1 , b = 0.12 and 7 = 0.131 (|McCulrough et, alJl200(i iHolman et"aH l20061. and 
find a\ = 95° and cry = 0.39 kms -1 for A = 0, and <j\ = 20° and cry = 1.8 kms -1 for A = 90°. Again, for this fairly 
central transit, the expected uncertainty depends strongly on A. 

For both TrES-1 and XO-lb the uncertainties in A and V are expected to be fairly large unless a substantial 
commitment of resources are expended to acquire many high-5/A spectra during transit, or the individual velocity 
uncertainties can be reduced substantially below the 15 ms -1 assumed here. This is because in both cases, the 
projected stellar rotation velocities are fairly small, and the transits are nearly central (small b). This gives rise to 
large uncertainties in V (for A ~ 0) or A (for A ~ 90°). 

In addition to these 5 systems, there are 5 other known tra n siting planets, all of which were originally identified as 
candidate s by the OGLE collabor ation iUdalski et al.l2002albllcll2003D. and subse quently con firmed with radial ve locity 
follow-up llKonacki et alJl2003albHKonacki et alJl2004 l2005UBouchv et al.ll200l iPont et alJl200l iMoutou et a"Tll200l 
iPont et al.H2005HBouchv et al.ll2005ajL The stars in these systems are all fainter than the five cases considered above 
(my ~ 15.5 — 17 versus my ~ 8 — 12), which will make the measurement of the RM effect much more challenging 
for these systems. Nevertheless, there is a good reason to expend this additional effort on at least some of the OGLE 
systems, namely, the "very hot Jupiters" (VHJs) with orbital periods less than 3 days. It has been suggested that 
these shorter-period planets discovered by OGLE form a disti nct population from the more numero us "hot Jupiters" 
with 3-4 d ay periods, based on differences in the ir typical mass l)Gaudi et alJl2005llMazeh et al.ll2005|) and their overall 
frequency l|Gaudi et aUl200fil iGould et, al.ll20f)l . This raises the interesting possibility that the VHJs arrived at their 
very close orbits through a different mechanism than the ho t Jupiters. For exam ple, the observation that the orbital 
distances of VHJs are nearly equal to twice their Roche radii l|Ford fc Rasiol20'06j) may indicate that these planets were 
emplaced v ia tidal capture and subsequent circularization. Such a scenario would b e the natural result of planet-planet 
scattering dRasio fc FordlH996l iFord et alJl2001l) . capture of free-floating planets l)Oaudill2003lk or Kozai oscillations 
l)Wu &: Murravl2003l D. Fabrvckv & S. Tremaine, priv. comm.). In the latter two scenarios, at least, one would expect 
these planets to generically have orbits with large misalignments with the spin axis of their parent stars. 

We can repeat our previous analysis to estimate the achievable uncertainty in A for the OGLE systems. Unfor- 
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Fig. 3. — The achievable total signal-to-noise ratio in RM waveform (left), and the signal-to-noise ratio in the the RM para meters A and 
V (right), as a function of b and A. Plotted is the "normalized" signal-to-noise ratio (i.e., after dividing by Qr; see Eg. 1181 for A (black 
lines) and for V (gray lines). The solid lines are based on our analytic formulae. The dotted and dashed lines are based on numerical 
computations with the more accurate formulas of Ohta et al. (2005); the dotted lines are for the case of no limb darkening, and the dashed 
lines are for linear limb darkening with e = 0.6. 



tunatcly, many of the parameters needed to provide an accurate estimate (particularly the impact parameter b and 
projected stellar rotation speed V) have not been measured, or have substantial uncertainties. We will therefore adopt 
approximate, fiducial values in order to provide an order-of-magnitude estimate, with the caution that the actual 
uncertainties could be significantly different. We will assume b = 0.5, 7 = 0.1 and e = 0.6. Only upper limits for the 
rotational velocities of the host stars have been reported, so we will adopt the largest allowed velocity of V = 5 kms -1 . 
Typical radial velocity precisions reported for these systems are a ~ 50 ms" 1 for ~ 45 minute exposures (using 
lOm-class telescopes), which we scale to ~ 100 ms -1 for 10-minute exposures. Assuming continuous observations 
during a typical transit duration of ^2 hr (i.e., twelve 10-minute exposures), we find Qr ~ 1.7 per transit, and so 
<j\ = 2Q i? 1 ~ 70°. With observations of 4 transits, one could measure A to within about 35°. While this is not nearly 
as good as the results that can be achieved for the systems with brighter host stars, even the finding that one of the 
VHJs has (say) A = —90° ± 35° would be quite interesting. 

4. DIFFERENTIAL ROTATION 

In the preceding calculations we neglected differential rotation, i.e., we assumed that the rotation speed of the star 
was independent of the stellar latitude. One might wonder whether or not the expected degree of differential rotation 
is detectable through transit observations. If so, it would be important both as a stellar astrophysics tool, and (if it is 
not modeled properly) as a possible source of ambiguity or covariance in the interpretation of RM data. 
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Differential rotation is an important phenomenon in stellar astrophysics because it is an observable manifestation of 
convective dynamics, and because it is intimately linked to the dynamo mechanism of stellar magnetic field generation 
and evolution. The Sun rotates ~20% faster at its equator than at high latitudes (see, e.g., Howard 1984 and references 
therein). Differential rotation has also been measured for other stars using star spots. In some cases, spots at different 
latitudes produce measurably different photometric periods (see, e.g., Henry et al. 1995, Rucinski et al. 2004, Herbst 
et al. 2006), and in other cases, the motion of star spots has been tracked via Doppler imaging (Collier Cameron et 



In principle, transiting exoplanets can be used instead of star spots. The advantages would be that planetary orbits 
are much more stable than star spot patterns, the transit chord can occur at any stellar latitude and can span a 
wide range of latitudes (unlike star spots), and the time scale of the transit (a few hours every few days) is more 
convenient than the month-long rotation periods of many stars. In addition, it would allow differential rotation to be 
measured on inactive and unspotted stars for the first time. A system of multiple transiting planets with different 
impact parameters has yet to be discovered, but would be a wonderful probe of differential rotation. But even with a 
single transiting planet, it is possible in principle to detect differential rotation through the RM effect if A ^ 0. The 
RM waveform effectively traces the projected surface speed of the star along a one-dimensional chord, and is sensitive 
to differential rotation if the chord spans a significant range of latitude. 

It is customary to parameterize observations of differential rotation according to a function such as V(i) = V cq (l — 

a sin 2 £), where V eq is the rotation speed at the stellar equator, I is the latitude on the surface of the star, and a is the 
differential rotation parameter (which is equal to the fractional difference in rotation speed between the pole and the 
equator). For the Sun, differential rotation is significant, a — 0.2. For stars with convective envelopes (T e // < 7500 K), 
the amount of differential rotation is thought to be correlated with V cq , such that more rapid rotators exhibit less 
differential rotation (Barnes 2005, Reiners 2006). For stars lacking exterior convection zones, the surface differential 
rotation is small and a ~ 0. 

If the surface of the star is in uniform rotation, then the surface speed at any point on the star depends only on the 
projected rotation speed of the star, V = Vs sin Is, and the distance from the the projected rotation axis, x p (Eq. |HJ. 
However, in the case of differential rotation, this symmetry is broken. The surface speed at a given point on the star 
also depends on the position parallel to the projected axis y p , as well as on the inclination of the rotation axis, Is- 
The apparent rotation speed at any point (x p ,y p ) of the stellar surface is 



Figure 0| shows the effect of differential rotation on the RM waveform. For this illustration, we used the parameters 
appropriate for HD 209458 (7 = 0.121, b = 0.5, and V cq = 5 krns" 1 ) and further assumed Is = 90°, so that the 
rotation axis is in the plane of the sky. We took the differential rotation parameter to be equal to the solar value of 
a = 0.2. Each column shows the results for a different value of A. The top row shows the waveform itself, and the 
bottom row shows the difference between the waveforms with a = and a = 0.2. 

For A = and Is = 90°, the effect of differential rotation is degenerate with changing the value of Veq, and hence 
differential rotation is not measurable. However, from Figure 0] and equation Q20JI. it is clear that for A 7^ or 
Is 7^ 90°, these two parameters are not degenerate, and in principle both parameters can be constrained by fitting 
a parameterized model to the velocity data. However, the requirement on the velocity precision is stringent. The 
amplitude of the deviations due to differential rotation generally scale with a, and are only a few m s _1 for a = 0.2. 

For a more quantitative analysis, we can apply the same procedure that we used in the previous section to estimate 
how well one can hope to estimate a using the RM effect. We estimate the uncertainties by nu merically computing the 
Fisher matrix (Eq. 115(1 using the analytic form for the RM effect given by lOhta et "al] l)2005|) , but accounting for the 
effect of differential rotation on projected surface velocity (Eq. |20J- We consider the uncertainties on four parameters 
a = (A, v cq , a, Is), and we assume that all other parameters have negligible uncertainties. The resulting uncertainties 
in the parameters a exhibit complicated relationships with A, b, and Is- We will summarize our results by focusing on 
the case a = 0.2, b = 0.5, 7 = 0.121, and e = 0.6 (approximating the HD 209458 system). Figure shows the expected 
uncertainties as a function of A for Is = 30°, 60°, 90°. We show the "normalized uncertainties," after dividing by Qr, 
where Qr is still given by equation l|18f) . but with a slightly altered definition of RM semiamplitude from equation 



First, we consider the case in which the stellar rotation axis is in the plane of the sky: Is = 90° (left column of Fig. 
|SJ). As seen in Figure 01 for this case and A = 0, a is degenerate with V eq , and it is not possible to determine these 
parameters independently. In fact, for Is — 90° and A = 0, neither of the differential rotation parameters (a, Is) are 
well constrained. We find that in general the uncertainty in A is degraded by covariances with the other parameters 
(particularly a). For example, we find a\ ~ SQr 1 f° r A ~ 0. This is ~ 50% larger than the uncertainty expected 
under the assumption of solid body rotation (Eq. I16f) . Thus without an a priori constraint on the allowed degree 
of differential rotation, measurements of A will be significantly compromised by an unknown amount of differential 
rotation. We explore the effects of a constraint on a in more detail below. 

From inspection of Figure we see that for Is — 90°, the most favorable value of A for measuring the overall system 
parameters (in the sense of minimizing the sum of the squares of the uncertainties of all four parameters) is A ~ 55°. 
For this value, we find <j\ = 4.1(5^ , av cq /V eq = 3.0(5^ , a a = 35(5^ , and <ji s = 52Q^ . Therefore measuring a to 
within 0.2 and Is to within 15° requires Qr ~ 200. Achieving this high a signal-to-noise ratio for a system similar to 



al. 2002). 
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©: K R = V cql 2 /(1 - 7 2 ). 
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Fig. 4. — The spectroscopic transit in the presence of differential rotation. The top row of panels shows the expected RM waveform for 
a transit in a system similar to HD 209458 (b = 0.5, 7 = 0.121, V cq = 5 kms ), for three different values of A. The stellar rotation axis 
is taken to be in the plane of the sky, and the differential rotation parameter was set to a = 0.2. The dotted lines show the corresponding 
waveforms for a = 0. The bottom row of panels show the difference between the a = 0.2 and a = cases. 



HD 209458b would require at least ~30 radial velocity measurements with uncertainties of ~ 2 ms 1 (although, of 
course, HD 209458 itself has proven to have a much smaller value of A). 

For Is = 60° (middle column of Fig. [5J the most favorable value is A w 70°. In this case, we find a\ = HQ^ 1 , 
<Jv eq /V eq = 28Qp , <J a — 31Qr i an d &i s — 42Q 7? 1 . The variances of the parameters are highly correlated. In 
particular, while the uncertainties in V eq and 1$ are individually high, they are correlated such that the uncertainty 
in V eq sin is considerably smaller. Measuring all of the parameters to better than 20% will require Qr > 200. 

We also computed the uncertainties for the same geometry as above, but for different values of Is and the impact 
parameter b. The right column in Figure shows the results for b — 0.5 and Is = 30°. Generally we find that the 
uncertainties are relatively large for Is < 30°, primarily because the amplitude of the RM effect is decreased by the 
sky projection. Grazing transits (b > 0.9) and central transits (b rj 0) provide poorer constraints, but for 0.2 < b < 0.7, 
the uncertainties are relatively insensitive to 6. 

We conclude that, for favorable geometries (30° < A < 70° and 0.2 < b < 0.7), it may be possible to detect 
differential rotation and to constrain a and Is to within ^20%, using very high signal-to-noise ratio observations of 
the RM effect (Q R > 10 2 ). 

On the other hand, for lower signal-to-noise ratio observations, we have shown that the covariances with the differen- 
tial rotation parameters (particularly a) will degrade the measurement of A. It is is therefore of interest to ask to what 
degree a modest a priori constraint on a will improve the expected uncertainties on A. Because such constraints are 
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Fig. 5. — The achievable signal-to-noise ratio in the RM and differential rotation parameters A, V eq ,a, Is, as a function of A for three 
values of Ig. We have assumed b = 0.5, 7 = 0.121, and e = 0.6, thus approximating the HD 209458 system. We have assumed differential 
rot ation parameter equal to the solar value of a = 0.2. Plotted is the "normalized" signal-to-noise ratio (i.e., after dividing by Qn; see 
Eg. 1181 . The solid lines show the case for no limb darkening, and the dashed lines are for linear limb darkening with e = 0.6. The dashed 
lines in the top row show our analytic formula for the uncertainty in A in the absence of differential rotation. The difference between 
the dashed and solid/dotted lines demonstrates the degradation of the uncertainty in A due to covariance with the differential rotation 
parameters. 



not easily incorporated into the Fisher matrix formalism, we instead determine the uncertainties using Monte Carlo 
simulations, as described in We first compute the expected uncertainties on the parameters with no constraints, 
and so verify our estimates based on the Fisher matrix formalism. 

Again adopting the parameters of the HD 209458b system (b = 0.52, 7 = 0.121, e = 0.6, A = -4?4, V eg sin/ S = 
4.70 kms' 1 , N — 14, and a = 4.1 ms _1 ), we find o\ = 2?7 for a — 0.2 (with a very weak dependence on Is). We then 
enforce a weak constraint on a by adding a penalty term to % 2 , of the form A% 2 = [(a — 0.2)/0.4] 2 . We find <j\ = 2?0, 
also with a weak dependence on Is- Thus, this mild a priori constraint reduces the expected uncertainty in A from a 
level that is 60% larger than the expectation under the assumption of solid-body rotation, to a level that is only 20% 
larger. At higher signal-to-noise ratios, the constraint on a has little effect. 

5. CONFIRMING TRANSITING PLANETS 

In this section we move from the regime of high S/N to the regime of low S/N, and investigate the utility of the 
RM effect in confirming the occurrence of transits in the most challenging cases. The photometric detection of transits 
offers one of the most promising ways to achieve the appealing goal of detecting terrestrial planets in the habitable 
zones of solar-type stars. That said, the challenges of robust detection of the photometric signal are daunting. A 
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transiting Earth-sized planet orbiting in the habitable zone of a solar-type star would reduce the stellar flux by only 
~1CP 4 , for a duration of only ~13 hr out of the year. Furthermore, the probability that a randomly oriented orbital 
plane happens to be close enough to edge-on to allow transits is only ~0.5%. The requirements for detecting these 
signals have convinced most researchers that it is necessary to conduct the search with space telescopes rather than 
ground-based telescopes. Several satellite experiments are slated for launch over the next few years that are desi gned 
to surv ey 10 4 — 10 5 m ain-sequence stars in order detect small transiting planets, including Corot l)Basrlinll2003[) and 
Kepler l)Borucki et al.ll2003j) . 
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Fig. 6. — Simulated spectroscopic signal of a transiting terrestrial planet in the habitable zone a solar-type star. The line-of-sight 
component of the stellar orbital velocity is visible as the sinusoid with a period of one year. The RM effect is visible as the spike at the 
origin, which is the time of transit. A close-up of the transit interval is shown in the lower panel. We have assumed b = 0, A = 0, e = 0.6, 
M = Mq, R = Rq, r = R<$, m = Af®, e = 0, P = 1 yr, and V = 5 kms" 1 . 

As with many such experiments, the most exciting discoveries from these surveys are likely to be those that are 
detected with the lowest signal-to- noise ratio (S/N). Furthermo re, the scaling of the number of detections as a function 
of the limiting (S/N) m i n is generally quite steep. For example. IGould et alJ l|2004j) demonstrated that the number of 
expected detections for Kepler scales as (S/N)~^? n . This implies that the number of detections is fairly sensitive to 
unanticipated degradations in (S/N) m - ln (from unexpected noise sources or other problems), and that the majority 
of the detections will occur very near the limiting (S/N) m i n . Since the limiting value of (S/N) m i n will generally be 
the minimum possible S/N for which detection is statistically possible, it will be difficult or impossible to use any 
further characteristics of the transit light curves themselves to distinguish true planets from false positives. If all of 
the information is required for mere detection, then there is not enough information for the accurate determination of 
multiple parameters. 
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Fig. 7. — Confirmation of Kepler planets. The symbols show the masses (or minimum masses) of known exoplanets versus their 
orbital semimajor axes. The black triangles are planets detected via spectroscopic orbits, red stars via transits, grey hexagons via pulsar 
timing, and purple circles via microlensing. Also shown are the Solar system planets. The green shaded region is where Kepler can detect 
(with S/N > 8) at least two transits by a planet orbiting a star with apparent magnitude ray = 12. The lines show the lower mass limits 
of planets detectable by Kepler that can be confirmed (with S/N > 5) via spectroscopic measurement of the stellar orbital over the course 
of a full period (red lines) and via measurements of the Rossiter-McLaughlin effect throughout a single transit (blue lines). We show these 
limits based on the capability of current instruments such as HARPS (solid lines), and the expected capability of a future 30m telescope 
(dashed lines). For this figure, we have assumed that the host star has the mass and radius of the Sun, that its projected rotation speed is 
5 kms -1 , and that all planets have the same mean density. Confirmation via the orbital velocity gives a lower mass limit in all cases, but 
it requires ~10 times more observing time, and greater instrumental stability (the spectroscopic orbit occurs over a full orbital period, as 
opposed to a single transit). 



Therefore, it will be essential to verify the low S/N candidates with additional observations. Unfortunately, since the 
photometric noise requirements are so stringent, it may prove very difficult to confirm these candidates with ground- 
based photometric observations. Space-based satellites such as HST or JWST l|Charbonneau et al.ll2004t iGould et alJ 
2004) may deliver the sensitivity and precision to confirm these detections, provided they are available. However, it 
would certainly be more expedient to have a good ground-based method to confirm candidates. 

The most desirable type of confirmation would be the spectroscopic detection of the stellar orbital velocity. This 
could be done from the ground, and the Doppler signal would also provide an estimate of the mass of the planet, and 
hence its mean density (when combined with the radius measurement from the photometric transit). The difficulty 
with the detection of the spectroscopic orbit is that the stellar reflex velocity would be quite small: the aforementioned 
Earth- like planet in the habitable zone of a solar- type star produces a ~9 cm s -1 wobble. Furthermore, the parent stars 
will be relatively faint (with an apparent magnitude my ~ 12 for Kepler), owing to the narrow-field, magnitude-limited 
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design of the survey observations. 

Can current setups detect the spectroscopic orbits of the systems with habitable planets that will be detected by 
Kepler? Before answering this question, we first estimate the parameters of the planets that can be detected with 
Kepler. The signal-to-noise ratio of a transiting planet with a semimajor axis a and radius r orbiting a star with a 
radius of R is 

£) T ^>- 1/2 (I)V^) 2 , pi) 

where TT is the total number of photons collected from a star of the appropriate apparent magnitude during the 
mission lifetime T. For the parameters of the Kepler mission, we find 

g) _„( « r 1/2 feV'Vi^V'V"*.^ ,, m 



NJ T vau/ \m®J V °- 75 

where we have assumed that R = Rq and that the planet has the same density as Earth: r = fl^m/M®) 1 / 3 . 

We now estimate the signal-to-noise ratio with which follow-up radial velocity measurements can confirm the Kepler 
detections. Assuming a circular orbit, and No radial velocity measurements are made with an uncertainty of a evenly 
sampled in orbital phase, then the total S/N of the orbital velocity signal is 

where we have defined 

Qo = VN^—. (24) 

(T 

The HARPS spectrograph ijPene et alJl2fl02h mo unted on th e 3.6m ESO telescope is representative of the state of 
the art in precision radial-velocity measurements. lLovis et alJ l)2005|) discuss observations of three dwarf stars with 
the HARPS setup. They take inventory of the contributions to the total error in HARPS measurements, including the 
photon noise, wavelength calibration, guiding errors, and stellar jitter (which in turn includes both stellar oscillations 
and granulation noise). They state that the uncertainties due to wavelength calibration and guiding errors can 
probably be im proved, and that it may be possible to average down the noise due to stellar jitter by using sufficiently 
long exposures l)Lovis i iTallEM but see lBouchv et all2005kl) . We w ill therefore focus on photon noise as the ultimate 
limiting noise source. Scaling from the results of lLovis et alJ l)2005l) . we can expect HARPS to achieve a ~ 1 ms" 1 
precision on a star with apparent magnitude my = 12 in a single 60 min exposure, depending on the spectral type of 
the host star. Assuming that No radial velocity measurements are made over the course of the year- long period of a 
terrestrial planet, this gives 

S ^ 0.6 V2 (JL) ^ (-£-) 1(r o,(^- 12 ) (D) , (25 ) 
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where we have assumed a primary mass of M = M©. Therefore, confirmation of habitable Earth-mass planets detected 
by Kepler will be difficult with current setups, even if stabilities of ~ 0.1 ms over a time scale of ~1 yr can be 
achieved. 

The requirements for the detection of the RM effect of a transiting habitable planet are generally less severe, both 
because the amplitude of the RM effect is larger than the orbital signal, as discussed in jj^J and because the time scale 
over which stability must be maintained is much shorter. Both of these points are illustrated in Fig. EI The S/N of 
the RM effect is approximately 



S 

~N ' = , - / " 



i(l -46 2 )cos 2 A + 6 2 



"I 1/2 

(26) 



Note that for A = 0, (S/N) R = Q R y/(l - & 2 )/3, whereas for A = 90°, (S/N) R = Q R b. Assuming a ~ 1 ms" 1 as 
above, M — M Q , R — i? Q , b — 0.5, and continuous measurements during a single transit (i.e. N R ~ 11.23 hr/1 hr), 
we find 

i) ^ 0, (r^r) (£=) ■ M 



NJ R ' VAU/ \ M ©/ V 5kms_1 / V 3 - 6m , 

Thus the S/N with which the RM effect is measured with only N R = 11 points during transit is roughly equivalent to 
the S/N with which the orbital Doppler shift is measured with No = 100 points. Furthermore, detection of the RM 
effect only requires stability at the ~0.4 ms" 1 level for ~11 hours, as opposed to a year. 

A downside to RM confirmation is that it must take place on the particular nights when a transit occurs, making 
the effort especially vulnerable to the vagaries of the weather. And of course, since the transit only represents a small 
fraction of the planet's orbit, it is possible to acquire many more data points on the orbital velocity curve than on the 
RM curve; it is possible to gain in the S/N by a factor of \J ira/R, where a is the planet's semimajor axis. However, in 
practice it will be difficult to fully realize this additional factor, because it is not possible to acquire data continuously 
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(due to telescope availability, weather, seasonal observability, and so forth), and because the stability requirements are 
harder to achieve for measurements over this longer time baseline. 

Nevertheless, the fact that the expected S/N is substantially less than unity even for continuous observations implies 
that it will only be possible to verify the most favorable candidates with current setups. Assuming that we can simply 
scale the photon noise uncertainties obtained by HARPS to larger apertures, we find that apertures of ~ 20 — 40 m 
will be required to confirm [with (S/N)r — 5] the habitable transiting planets orbiting my =12 primaries detected 
by Kepler. 

Figure [7| summarizes the prospects for the detection of planets via Kepler and the prospects for the confirmation 
of those planets via spectroscopic detection of either the orbital velocity or the RM effect. We show the region 
of parameter space in the m — a plane in which Kepler can detect transiting planets orbiting stars with apparent 
magnitude my — 12, assuming at least two transits are required for detection, and (S/N)t > 8. By design, Kepler 
is (just) sensitive to Earth-mass planets in the habitable zone. We also show the region of parameter space in which 
planets can be confirmed (with S/N > 5) by detecting either the orbital velocity or the RM effect using current 
facilities, and using the expected capabilities of future 30m telescopes. 

6. SUMMARY AND DISCUSSION 

The transit of an exoplanet is a rare and information-rich event that should be exploited in every possible way. In 
this paper we have examined one aspect of transit physics, the Rossiter-McLaughlin effect, whose origin is the rotation 
of the parent star. Our goal has been to assess the ability of near-term observations to exploit the information in the 
RM signal. 

We have calculated the achievable accuracy with which one can measure the key parameter A, with which one can 
assess the degree of alignment between the stellar spin axis and the planetary orbit, and have assessed the prospects for 
the currently known sample of transiting exoplanets. One might ask further: how accurately must A be known in order 
to enable interesting theoretical advances, such as the discrimination among different planetary migration theories? 
In other words, how good is good enough? This question does not seem to have been addressed in the literature on 
planet formation theory, and is an appealing topic for future research. As suggested in § 3, for a system with a very 
hot Jupiter, even a crude measurement with o\ = 30° would be of interest. As a more general benchmark, one would 
like to achieve at least enough accuracy to tell whether a given system is similar to the Solar system, or not. To date, 
the only system that has been observed with enough precision for this task is HD 209458, and the result is that it is 
Solar- like with a small but significant misalignment of A « —4° (Winn et al. 2005). A system could be different cither 
by having a larger misalignment (A ^ 10°), for which an accuracy of a few degrees or more would suffice, or by having 
a more perfect alignment, for which an accuracy of ~ 0.1° would be needed. 

We have also explained and quantified the potentially important role of the RM effect in the confirmation of transiting 
planet candidates. It is a difficult task, but in many cases it will be easier than the task of detecting the orbital velocity 
of the parent star. The comparative advantages of RM confirmation are that the velocity amplitude is generally larger, 
and the acceleration is certainly larger. The full range of RM variations occur during a single transit (typically a 
single night or two) whereas the orbital velocity variations occur over the full orbital period (a year, for the interesting 
case of an Earth analog). A disadvantage of RM confirmation is that one does not learn the planetary mass; it is not 
a dynamical measurement, per se, but rather an alternate method of verifying that a portion of the stellar surface 
is periodically occulted. Furthermore, RM observations must take place during transits, whereas the spectroscopic 
orbital data can be obtained on any other nights. In addition, the ease or difficulty of the RM method depends on the 
rotation rate of the star, which will vary from system to system. 

Another interesting question for future research is whether or not the RM effect could ever profitably become 
the basis of a planet search technique, rather than simply a follow-up technique. Ohta et al. (2005) suggested that 
transiting planets might announce themselves through large outliers in data bases of radial velocity measurements. The 
relevant efficiency factors, selection effects, and detection algorithms for such an undertaking have yet to be worked 
out. Even less certain are the prospects for an RM-based transit survey, in which stars are spectroscopically monitored 
specifically for transits. The obvious drawback is that typical high-resolution spectrographs examine only one star at 
a time, whereas photometric transit surveys observe thousands or even millions of stars at a time. There are plans for 
multiplexed Doppler spectrographs in the near future (Ge et al. 2004), but the multiplexing factor is only ~50. 

On the other hand, the RM effect is perhaps uniquely sensitive to small planets around hot and rapidly rotating 
stars. The amplitude of the RM effect depends on the stellar radius and projected rotation speed, according to Eq. @, 
but these factors are not completely independent. Among main-sequence stars, both the stellar radius and the typical 
rotation speed increase with stellar mass (or temperature), with an especially dramatic rise near the F5 boundary. The 
typical rotation speed of GO stars is ~10 km s _1 , and for F0 stars it is ^100 km s _1 (see, e.g., Tassoul 2000). This is 
related to the strong variation in the depth of the convective zone across this same boundary. We have attempted to 
illustrate the net effect of these astrophysical correlations on the strength of the RM signal, in the following manner. 
We considered main-sequence stars ranging in spectral type from 05 to GO, and imagined that a Jupiter-sized planet 
makes an equatorial transit of each star. 4 We calculated Kr using estimates of the stellar radius from Cox et al. (2000), 
and estimates of typical rotation speeds from Tassoul (2000). The results are shown in Fig. [S] The maximum signal 
is achieved for late A stars. Of course, the achievable velocity accuracy is also a function of stellar type, as it depends 

4 Later-type stars are not amenable to even this crude analysis, as there is a very wide dispersion in velocity (and achievable velocity 
accuracy) that is linked to age and stellar activity through the relations of Skumanich (1972) and Noyes et al. (1984). 
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on the number and width of absorption lines available for velocity determination. We have not attempted to quantify 
the velocity accuracy because in this regime it would certainly be advantageous to abandon the description of the RM 
effect as a net Doppler shift and model the distortion directly, as is done in Doppler-imaging algorithms. 



700 P 




Spectral type 



Fig. 8. — Rough estimate of the amplitude of the RM effect of a transiting Jovian planet as a function of the spectral type of the parent 
star. The effect reaches a maximum for late A stars, which is a compromise between the faster rotation rates and the larger radii of hot 
stars. 

As with so many other areas of stellar and planetary physics, interest in the Rossiter-McLaughlin effect has been 
re-invigorated by the discovery of exoplanets. It has been nearly 100 years since the first observations of this effect, 
and yet we believe that there are still important applications waiting to be developed. 



We are grateful to Thomas Beatty, Tim Brown, Bill Herbst, Matt Holman, Ed Turner, and Yasushi Suto for helpful 
discussions. B.S.G. was supported by a Menzel Fellowship from the Harvard College Observatory. 
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